Co-doped Ceria: Tendency towards ferromagnetism driven by oxygen vacancies 
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We perform an electronic structure study for cerium oxide homogeneously-doped with cobalt im- 
purities, focusing on the role played by oxygen vacancies and structural relaxation. By means of 
full-potential ab-initio methods, we explore the possibility of ferromagnetism as observed in recent 
experiments. Our results indicate that oxygen vacancies seem to be crucial for the appearance of 
a ferromagnetic alignment among Co impurities, obtaining an increasing tendency towards ferro- 
magnetism with growing vacancy concentration. The estimated couplings cannot explain though, 
the experimentally observed room-temperature ferromagnetism. In this systematic study, we draw 
relevant conclusions regarding the location of the oxygen vacancies and the magnetic couplings in- 
volved. In particular, we find that oxygen vacancies tend to nucleate in the neighborhood of Co 
impurities and we get a remarkably strong ferromagnetic coupling between Co atoms and the Ce^"^ 
neighboring ions. The calculated magnetic moments per cell depend on the degree of reduction 
which could explain the wide spread in the magnetization values observed in the experiments. 



I. INTRODUCTION 

The terrain of the so called Diluted Magnetic Oxides 
(DM0) is currently being explored with a strong drive to 
find Room Temperature (RT) ferromagnetic (FM) mate- 
rials for possible technological applications^. DM0 are 
experimentally obtained by doping the oxide matrix with 
a small amount of magnetic Transition Metal (TM) ions. 
This procedure may introduce ferromagnetism in other- 
wise non magnetic materials and opens up the feasibility 
of applications that range from spintronics to magneto- 
optical devices. 

Since the pioneering works of Dietl et al.^ in Mn-doped 
ZnO, RT ferromagnetism was observed in various doped 
oxide hosts- with different TM ions and doping concen- 
trations. Some DM0 are FM insulators while others are 
semiconductors and, so far, there is not an obvious link 
between the ferromagnetic behavior and the conduction 
properties. The samples are sensitive both to prepara- 
tion methods and growth conditions and this observation 
supports the idea that ferromagnetism might be closely 
linked to defects and/or to the presence of oxygen vacan- 
cies. Interesting for spintronic applications are critical 
temperatures being far above room temperature. 

Among DM0, Co-doped Ce02 (Coa;Cei_a;02) has 
attracted particular interest due to ferromagnetic be- 
havior observed well above RT for low Co-doping 
concentration^i^iii^. Moreover, ceria is a transparent 
and high dielectric constant rare-earth oxide, whose flu- 
orite structure matches well with that of silicon. It 
keeps the crystal structure, both, under doping and upon 
the formation of oxygen vacancies, thus promising a 
good integrability for spintronic devices even for non- 
stoichiometric compounds. However, there is controversy 
among the experimental results regarding critical tem- 



peratures and the values of the magnetic moments, as a 
function of doping concentration. As a matter of fact, a 
wide range of magnetic moments has been reported^s^i^ 
and this has raised concerns about the intrinsic nature of 
the FM properties of these materials due, for example, 
to Co secondary phases, heterogeneities or even contam- 
ination. 

The presence of oxygen vacancies has been considered 
as a possible factor affecting the FM response in the case 
of films^ but it is not yet clear whether it induces Co clus- 
tering and/or promotes magnetic ordering. It has been 
recently argued that even without introducing magnetic 
impurities, the presence of oxygen vacancies in Ce02 
could stabilize ferromagnetisn*^. 

There is experimental evidence available, obtained us- 
ing several technique a^^i^^'^^i"'^^'^"^ , indicating that the 
charge left behind by the oxygen vacancies gets localized 
near some Ce atoms driving the Ce^"*" to Ce'^+ reduction. 
In the case of undoped Ceria surfaces, this localization 
has been recently addressed by Ganduglia-Pirovano et 
al}-^. The charge localization happens both in doped as 
well as in undoped Ce02-s, and has to be properly taken 
into account since it might affect atomic relaxation and, 
in turn, the magnetic behavior. It is important to remark 
that in order to perform a theoretical study of the mag- 
netic properties of Co-doped reduced ceria, it is critical 
to relax the crystal structure. 

Previous ab-initio DFT+U calculations have reported 
a detailed description of the relaxation processes in re- 
duced ceria upon doping with transition metal atoms 
such as Zr— and Pdii. It results that relaxation is dif- 
ferent in Zr-doped systems as compared to the Pd-doped 
ones, because the electronic distribution, which changes 
upon reductionii, depends on the size and the chemical 
characteristics of the dopant. 
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Our goal in this paper is to explore the role played by 
oxygen vacancies on the magnetic properties of a bulk 
ceria matrix as a function of Co-doping concentration. 
We perform Local Density Approximation (LSDA)+U 
calculations^^ considering bulk unit cells with substitu- 
tional replacement of Co atoms into cerium sites"'^^. In 
this work we take into account, in the relaxation process, 
the presence of Ce'^"'' ions induced by reduction. We also 
consider the evolution of the magnetic interactions with 
the number of oxygen vacancies per Co atom, which (to 
the best of our knowledge) has not yet been reported for 
doped ceria. 

We organize the paper as follows: In Sec. HIl and Sec. 
mil the computational method is described. In Sec. IIVI 
and Sec. |Vl the effect of introducing oxygen vacancies in 
Co-doped ceria is studied. In Sec. IVI[ the magnetic and 
electronic properties as a function of the dopant concen- 
tration and oxygen deficiency are explained and finally, 
in Sec. IVIH the results of this work are discussed. 



II. METHOD OF CALCULATION 

Density functional theory (DFT)^ calculations are 
performed within the Local Density Approximation 
(LDA)^J. We use the full potential augmented plane 
waves method as implemented in the Wien2k code^^, 
where the space is divided into muffin tin (MT) spheres 
around the atoms and an interstitial region. Plane waves 
are used to describe the region outside the spheres. The 
mufhn tin radii used are R^fj^ =1.9 au, R^fj^ =1.6 au, 
^MT =2.3 au. The number of plane waves in the intersti- 
tial region is set by the cut-off parameter IlKmax=7 (or 
6 for the biggest cells) where R is the minimum Rmt in 
the corresponding cell. To attain the desired convergence 
precision, we use a 7 X 7 X 7 k-mesh in the Brillouin zone 
for a Co-doping concentration of x=12.5 % and a 5 X 5 
X 5 mesh for x=6.25 %. 

In the cases where we have charge localization on the 
4/ states of Ce (which we call Ce'^"'"), the structures 
are first allowed to relax taking these 4/ electrons as 
core levels^^, until forces on atoms are below 1 mRy/au. 
Then, the electronic structure of the fully relaxed super- 
cell is obtained using the LSDA-f U approach switching 
on the local Coulomb interaction in the Ce'^^ ions^"^ with 
Ue//=U - J = 6 eV— . The relaxation processes involved 
will be described in detail in Sec. |Vl 



III. SUPERCELL CALCULATIONS 

The crystal structure of Ce02 is fluoritc with an exper- 
imental lattice parameter of 5.411 A. Each Ce^+ cation is 
coordinated to eight nearest neighbors and in turn, 
each is tetrahedrally coordinated to Ce*+ cations. 

Previous theoretical results for doped ceria indicate 
very small changes in the lattice constants upon doping 



with different transition metals, while preserving the cu- 
bic symmetrji^'^^. In particular, when doping with Zr, a 
0.2% contraction in the lattice parameters is obtained^ 
and a 0.04% reduction, when doping with Pdr^. Tak- 
ing this fact into account, we consider the experimental 
lattice constants of unreduced Ce02 for all the systems 
under study. However, as we are interested in under- 
standing the effect of both, magnetic impurities and oxy- 
gen vacancies, it is essential to relax the internal atomic 
positions, as mentioned in the previous section. 

For a concentration of Co given by x=6.25 %, we 
use a 2x2x2 BCC-type supercell, and for x=12.5 %, a 
2x2x2 FCC-type one. The supercells are built out of the 
conventional 12-atoms cubic unit cell of Ce02 and are 
schematically shown in Fig. [T] For the unreduced com- 
pounds these supercells contain 48 and 24 atoms, respec- 
tively. The nearest neighbor distance between Co atoms 
is 9.37A in the BCC cell and 7.65 A in the FCC one, so 
that these atoms can be safely considered as impurities. 

For each concentration we study configurations with 
zero, one or two vacancies per Co impurity atom. When 
an oxygen vacancy is created in pure bulk Ce02, two 
electrons are left in the system and the vacancy is sur- 
rounded by four next nearest neighbor Ce atoms. As 
mentioned before, several experimental results seem to 
indicate that these extra electrons occupy localized 4/ 
states on two Ce sites. There are then, two types of Ce 
atoms: Ce^+ ions in a 4/°-like configuration and Ce'^+ 
ions in a Af-^ one, that need to be considered theoretically 
(see for instance Ref. |2^ andi26i, and references therein). 

When Co impurities are introduced, the system gains 
in complexity. X-ray experiments suggest that in oxy- 
gen deficient samples the valence of Co is 2-1- (see Refs. 
Mlslfl^). Regarding Ce, the presence of Ce'^"'" has been 
reported in Refs. ItHT^. In particular, in Ref. [13, for a 
sample with a 4.5 % Co concentration, the results seem 
to indicate that around 10% of the Ce atoms are Ce"^"*" 
while the others are Ce^+. 

Taking the above into account, in this work we con- 
sider a variety of situations, depending on the number 
of oxygen vacancies per impurity Co atom, as it will be 
detailed in the following section. 



IV. LOCATION OF THE OXYGEN VACANCIES 

The first question which arises when introducing va- 
cancies into the Co-Ce02 system relates to their location 
within the ceria matrix: what is energetically more fa- 
vorable, to have oxygen vacancies close to or far from 
the Co impurities? To answer this question we first place 
one vacancy in a site nearest neighbor to a Co atom and 
relax the structure for the case x=12.5 %. We then put 
the oxygen vacancy far away from the Co atom and af- 
ter relaxing the structure, compare both results. In view 
of the experimental evidence mentioned in the previous 
section and the model suggested by Vodungbo et al^, we 
consider that the extra electrons left behind by the va- 
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FIG. 1: (Color on-line) SuperCells considered for two diflterent 
Co-doping concentrations. Top: 2x2x2 BCC type supercell in 
which one of the 16 cerium sites is replaced by a Co atom. 
This supercell contains 48 atoms. Bottom: 2x2x2 FCC su- 
percell where one of the 8 cerium sites is replaced by a Co 
atom. This is a 24 atom-supercell. 



cancy, migrate to the Co atoms leaving all the Ce ions as 
Ce^+, for both vacancy locations. The difference in en- 
ergy between the two considered vacancy configurations 
is 0.55 eV per Co atom in favor of the vacancy being close 
to the impurity (see Fig [21(b). 

We also consider a larger concentration of vacancies, 
namely, two per Co impurity atom and perform the cal- 
culations for x=6.25 %. We take into account two possi- 
bilities: a configuration where the two oxygen vacancies 
are close to the Co atom as suggested in Ref. d (see 
scheme in Fig. [11(c)) and a configuration where one oxy- 
gen vacancy is close to the impurity while the other is far 
away from it. Comparing the total energies of both con- 
figurations, we find that the two oxygen vacancies prefer 
to be close to the Co atom by 0.62 eV per Co atom. 



It is clear from these results that cobalt has a strong 
tendency to nucleate oxygen vacancies in ceria. In the 
following sections we discuss the relaxation process, the 
electronic structure and the magnetic couplings for the 
most favorable spacial distribution of vacancies. 
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FIG. 2: (Color on line) Schemes of the lowest energy con- 
figurations of oxygen vacancies within the ceria matrix for: 
(a) Unreduced system, (b) One vacancy per Co atom, (c) 
Two vacancies per Co. The white square indicates an oxygen 
vacancy. 



V. IONIC RELAXATION 

A. Zero oxygen vacancies 

As discussed in the previous section, for the doped sys- 
tem without oxygen vacancies, both Co and Ce are ex- 
pected to be -1-4 and therefore we use conventional LSDA 
calculations to relax the supercells for the two impurity 
concentrations considered. The obtained relaxations are 
qualitatively similar: The eight oxygen atoms which sur- 
round the impurity get symmetrically closer to it by 0.09 
A (for x=6.25%) and by 0.07 A (for x=12.5%). Re- 
garding the nearest cerium atoms, they move towards 
the impurity by around 0.04 A for x=6.25%, while they 
do not relax due to symmetry reasons for the largest con- 
centration considered. 



B. One oxygen vacancy per Co atom 

The valence configuration is depicted in Fig. [H(b). 
As there are no localized 4/ electrons in the cell for this 
vacancy concentration, we use again conventional LSDA 
calculations and focus only on x=12.5%. 

After relaxation, the seven oxygen atoms, which stay 
around the impurity after creating the vacancy, undergo 
a displacement of around 0.06 A in average: one oxygen 
moves towards Co considerably (~0.3 A), three oxygen 
atoms get closer to it by ~0.1lA and the last three 
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ones move outwards by '-^0.07 A. Regarding the nearest 
neighbors of the vacancy, the nearest three Ce ions and 
the Co impurity move outwards with respect to it (the 
Ce atoms displace by 0.1 A and the Co atom by 0.28 
A) while the six nearest neighbor oxygen atoms move 
towards the vacancy by 0.15 A on the average. 



Zr-doping, the charge left behind by the vacancy local- 
izes on two neighboring Ce atoms, while for Pd and Co 
it localizes both at the impurity site and at its nearest 
neighbor oxygen atoms. When doping with Pd, the dif- 
ference in relaxation might be due to the different size of 
the dopant atomic radii. 



C. Two oxygen vacancies per Co atom 

We consider the most stable vacancy location, namely, 
the two vacancies being near the impurity atom. As men- 
tioned above, in this case, the extra electronic charge left 
behind by the vacancy, localizes at two Ce atoms, which 
are nearest neighbors to Co. 

After relaxation, the six oxygen atoms surrounding the 
transition metal impurity move towards it. For x=6.25 
%, two of them displace by ~0.27 A while the other four 
do it by ~0.15 A. The oxygen atoms which are nearest 
neighbors to the Ce'^"'" ions move away from them in av- 
erage about 0.1 A. This expansion of all oxygen atoms 
around the Ce'^+ has to do with the coulomb repulsion 
due to the extra charge localized at that Ce site. In con- 
trast, this is not observed around the Ce*"*" ions, where 
some of the nearest neighbor oxygen atoms get closer to 
these ions while others move away, but without under- 
going a net expansion of the oxygen-cloud. Qualitatively 
similar atomic relaxations are obtained for x=12.5%. It 
is clear that the introduction of vacancies produces an 
important rearrangement of the structure. 

D. Comparison with other dopants 

It is interesting to compare the results obtained when 
doping with Co with those appearing in the litera- 
ture for other dopants such as Zr or Pd. The avail- 
able data include cases considering just one vacancy per 
impurityi^iii. 

For unreduced ceria, relaxation in the presence of 
Zr impurities is similar to that of Co: the 8 nearest 
neighbor oxygen atoms and the three nearest Ce^^ ions 
move towards the impurity by ~0.1 A and ^--^0.036 A , 
respectivelyi^. This is due to the fact that Zr and Co 
have smaller atomic radii than the Ce^+ ion. When dop- 
ing with Pd, which has an atomic radius similar to the 
ionic radius of Ce^+, the displacements towards the im- 
purity are smaller than in the case of Zr and Coii. 

In the presence of one oxygen vacancy, relaxation for 
Co-doping resembles neither the Pd nor the Zr dopant 
cases. Co and Zr move away from the vacancy while Pd 
moves towards it. All nearest neighbor oxygen atoms get 
closer to the vacancy when the dopants are Co and Pd. 
In the case of Zr, most of those oxygen ions move towards 
the vacancy except for one oxygen that moves away from 

The above mentioned differences in the relaxation pro- 
cesses might be due to different reasons: In the case of 



VI. MAGNETIC AND ELECTRONIC 
PROPERTIES OF CERIA AS A FUNCTION OF 
CO DOPING AND REDUCTION. 

A. Magnetic coupling between cobalt and Ce'^'^ 
ions 

Before going into the analysis of the evolution of the 
electronic and magnetic properties of ceria as a function 
of dopant and vacancy concentration, it is interesting to 
determine the nature of the magnetic coupling between 
Ce^+ and Co. From the above discussion, it is clear that 
this coupling takes place when two vacancies per Co are 
present. To get insight into this problem, we consider the 
case of x=12.5% to compare the energies of FM and AFM 
spin- alignments of Ce^+ with respect to Co, as depicted 
in Fig. [31 The energy difference between these config- 
urations is 120 mcV per Co atom in favor of the FM 
alignment. Cobalt and Ce'^^ ions show, then, a strong 
ferromagnetic coupling. 
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FIG. 3: (Color on line) Scheme of the two configurations con- 
sidered for calculating the magnetic coupling between Co and 
Ce^+ ions (Top, FM. Bottom, AFM). The FM configuration 
is favoured by an energy difference of 120meV per Co atom. 



B. Electronic structure and magnetic coupling 
among dopants 

Pure Ce02, without impurities or oxygen vacancies, 
is a non-magnetic insulator. The occupied bands are 
mainly of O 2p character and the unoccupied ones are 
mostly of Ce 4f character (See Fig. d]). When impurities 
are introduced in the unreduced system, defect states ap- 
pear due to majority and minority contributions coming 
from the d states of the substitutional Co atoms. For 
both impurity concentrations, and no oxygen vacancies. 
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the magnetic coupling among Co ions is slightly antifer- 
romagnetic. In Table HI we show the energy differences 
among FM and AFM configurations, A E =Efm -^af 
for both Co-concentrations and different number of va- 
cancies per Co. 



case 


AE (x=6.25%) 


AE (x=12.5%) 


vac 


+0.4 


+2.5 


1 vac 


0.0 


-0.5 


2 vac 


-3.1 


-12.5 



TABLE I: AE —EpM - Eaf, energy difference per Co ion (in 
meV) between FM and AF solutions, to estimate the magnetic 
coupling between Co atoms, for different cobalt concentration 
X and number of vacancies. 

When oxygen vacancies are introduced, a tendency 
towards FM Co-Co coupling appears. The energy dif- 
ference A E changes sign thus favoring a growing FM 
alignment, which increases with, both, vacancy and Co 
concentration. For the largest Co concentration, the 
vacancy-induced ferromagnetic Co-Co coupling, is clearly 
established. 

Experimental evidence seems to suggest that there are 
no metallic Co clusters or secondary phase formation in 
Co-doped ceriEfc^. In this work, assuming then, that 
Co impurities are diluted and homogeneously distributed 
within the ceria matrix, a tendency towards ferromag- 
netism is obtained. This tendency increases with oxygen 
deficiency. However, the energy differences obtained can- 
not explain the high Curie temperatures observed in the 
experiments. Thus, the underlying mechanism for the 
observed room temperature ferromagnetism, remains un- 
clear. 

To compare the effect that the introduction of vacan- 
cies has on the electronic structure of doped ceria, we 
show the densities of states corresponding to the FM solu- 
tions. As it can be observed in Fig. [51(a) and Fig. [61(a), 
for the two unreduced Co concentrations considered here, 
two broad impurity peaks develop at the bottom of the 
majority valence band of 0-p character. The most im- 
portant contribution to the peak appearing at the top of 
the majority valence band, comes from the oxygen atoms 
which are nearest neighbors of the Co impurities. The or- 
bitals of these oxygen atoms are magnetically polarized 
due to hybridization with the Co d impurity orbitals. 
Two localized minority impurity peaks lie in the gap of 
ceria. For both Co concentrations, the Fermi level falls 
within the first minority peak, which lies at the top of the 
oxygen minority p band. The second and unoccupied mi- 
nority impurity peak appears well inside the original gap 
of ceria, far from the band edges. 

In order to better follow the evolution of the impu- 
rity peaks with the introduction of vacancies, we show 
in Fig. [71 the projected local densities of states on Co 
for x=6.25%. The trends are similar for x=12.5 %, with 
more delocalized impurity bands due to larger Co-Co hy- 
bridization. 



When oxygen vacancies are introduced, Co''"'" reduces 
to Co2+. In Fig. [51(b) and (c) and Fig. [i(b) and (c), the 
corresponding total densities of states and the Co local 
densities of states are shown. Due to the Coulomb re- 
pulsion the large spectral weight of Co d character, orig- 
inally lying at the botton of the majority valence band 
goes into two new majority peaks, which can be observed 
above this band. There are now three majority Co im- 
purity peaks lying below the Fermi level and above the 
valence band in these reduced configurations for the two 
Co concentrations considered. 

When two oxygen vacancies per Co are present, the 
extra charge, which localizes on two Ce^+ sites, gives 
rise to the occupied majority peak of 4f character, which 
can be observed in Fig. [51(c) and Fig. ^{c). 
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FIG. 4: (Color on line) Total density of states for pure Ce02 
(black line) and partial one for Ce (orange). 



C. Evolution of the magnetic moments 

In the non reduced impurity systems, the total mag- 
netic moment of the FM cells is mainly due to the po- 
larization of Co and its nearest neighbor oxygen atoms. 
The AFM solution shows a distribution of the magnetic 
polarization among the atoms which is not different from 
the one obtained for the FM solution. The hybridization 
effects are, namely, very local. The eight nearest neighbor 
oxygen atoms of Co become ferromagnetically polarized 
with respect to the impurity, for both Co Concentrations. 
See Tables [III and [ml 

When we introduce one vacancy, the two extra elec- 
trons which are left behind give rise to a decrease in the 
total magnetic moment of the ferromagnetic cell. The to- 
tal moment goes from 4.76 /is (for x=6.25%) and 4.60/is 
(for x=12.5%) to 3 /is. This decrease in the value of the 
magnetic moment per Co atom of almost 2 fiB points to- 
wards an effective Co valence which goes from +4 to +2. 
The extra charge left behind by the oxygen vacancy local- 
izes, then, in the electronic cloud surrounding the transi- 
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FIG. 5: (Color on line) Density of states for Cei-xCox02~s 
with x=6.25% for (a) unreduced system (5=0), (b) one va- 
cancy per Co (5=0.0313), and (c) two vacancies per Co 
(5=0.0625). Total DOS (black line), partial DOS for cobalt 
(blue) and cerium (orange). In (a) and (b), there is no charge 
localization on any Ce site. In (c) the partial DOS of a Ce'^"'" 
is shown. 



FIG. 6; (Color on line) Density of states for Cei-xCox02~s 
with x=12.5% for (a) unreduced system (5=0), (b) one va- 
cancy per Co (5=0.0625), and (c) two vacancies per Co 
(5=0.125). Total DOS (black line), partial DOS for cobalt 
(blue) and cerium (orange). In (a) and (b), there is no charge 
localization on any Ce site. In (c) the partial DOS of a Ce^+ 
is shown. 



tion metal ion, which is covalently shared by this ion and 
the seven remaining nearest neighbor oxygen atoms. 

Finally, when there are two vacancies per Co, the mag- 
netic moment per magnetic ion is 5 /is, for both Co con- 
centrations. From the four extra electrons per Co atom 
left behind by the vacancies, two contribute to decrease 
the magnetic moment of the cloud formed by Co and its 
neighboring oxygen atoms with respect to the unreduced 
situation. The other two electrons get localized on two 
Ce^+ sites close to the vacancies, giving rise to the in- 
crease in 2 fj,B of the magnetic moment of the cell per Co 
atom which goes from 3 fiB to bfiB- See Tables HIl and 

imi 

The values of the total magnetic moments of the fer- 
romagnetic cells also reveal that the introduction of va- 
cancies changes the band character of the doped systems. 
The unreduced impurity systems show a non integer cell 
magnetic moment, characteristic of a metallic-like band 
structure, while the introduction of vacancies gives rise 
to half-metallic-like features, as for instance the integer 
value of the total cell spin moments. 

It is interesting to remark, again, that the total mag- 
netic moment per cell strongly depends on the degree of 
reduction. This might explain the wide spread in the 
magnetization values obtained in the experiments as re- 



ported in the literature^i^i^i^. 



case 


/IT 


f^Co 


nn 

Mo 


^i{Co + o"") 


/i(Ce') 


vac 


4.76 


2.94 


1.28 


4.22 


0.04 


1 vac 


3.00 


2.50 


0.42 


2.92 


-0.17 


2 vac 


5.00 


2.49 


0.37 


2.86 


0.95 



TABLE II: Magnetic moments (in fis) for x=6.25 %: iiq" 
means the sum of the magnetic moments of all oxygen atoms 
nearest neighbor to cobalt. fi{Co + O"") means the mag- 
netic moment of the cloud formed by the Co and its nearest 
neighbors. fi{Ce') is the one corresponding to a cerium atom 
close to the vacancy. For the 2 vacancies case, there are two 
reduced Ce per cobalt. 



VII. DISCUSSION AND CONCLUSIONS 

In this work we have studied the effect of oxygen va- 
cancies on the magnetic and electronic properties of Co- 
doped ceria. We have shown that the introduction of 
oxygen vacancies is essential to drive a ferromagnetic cou- 
pling among the dopant impurities. Actually, in the ab- 
sence of vacancies there is a slight tendency towards an 
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FIG. 7: Local density of states on the Co impurity site, with 
increasing oxygen vacancy concentration for Co concentration 
x=6.25%. (a) without vacancies, (b) one vacancy per Co, (c) 
two vacancies per Co. 



case 


fJ,T 




nn 

Mo 






vac 


4.60 


2.90 


1.28 


4.18 


0.04 


1 vac 


3.00 


2.50 


0.40 


2.90 


-0.03 


2 vac 


5.00 


2.41 


0.40 


2.81 


0.99 



TABLE IIL Same caption as Table HH for x=12.5 %. 



antiferromagnetic coupling that changes towards a ferro- 
magnetic one, when the vacancy concentration increases. 

For two vacancies per Co, part of the left-behind charge 
localizes in Ce ions sitting in the vicinity of the impurity, 



turning them into magnetic Ce^+ ions. We find that the 
magnetic coupling among the Co impurities and the Ce'^"'" 
ions, is strongly ferromagnetic. This coupling provides a 
large magnetic moment, localized close to the impurity. 

Even if Co impurities in unreduced ceria show large 
local magnetic moments, the interaction among impuri- 
ties is small for the two dopant concentrations studied. 
The same can be said about the systems with only one 
vacancy per impurity atom. The ferromagnetic coupling 
among impurities in the presence of two vacancies per Co 
increases when going from an impurity concentration of 
6.25 % to 12.5 %. 

We conclude that, although there is a clear tendency 
towards ferromagnetism in Co doped ceria when oxygen 
vacancies are present, the obtained FM couplings cannot 
explain the high Curie temperatures observed in the ex- 
periments. On the other hand, the large values of the 
magnetic moments per Co atom observed can be under- 
stood if the Co impurities, the vacancies and the nearest 
Ce'^+ ions build complexes of the type considered in this 
work. To understand the simultaneous appearance of 
the two phenomena, namely Curie temperatures above 
RT and a wide spread in the large values of the magnetic 
moments per Co dopant, it might be necessary to con- 
sider eventual inhomogeneities and/or to introduce other 
magnetic interaction mechanisms which are not present 
in the exchange interactions taken into account within 
the framework of these DFT calculations. 
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